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ABSTRACT 

We investigate the variability properties of blazars in the GeV band using the data of the Fermi/LAT 
telescope. We find that blazars exhibit variability on the scales down to the minimal timescale 
resolvable by Fermi, which is a function of the peak photon count rate in the LAT. This implies that 
the real minimal variability timescales for the majority of blazars are typically shorter than those 
resolvable by the LAT. We find that for several blazars these minimal variability timescales reach 
those associated to the blazar central engine, the supermassive black hole. At the same time, none of 
the blazars exhibits variability on the timescale shorter than the black hole horizon light crossing time 
and/or the period of rotation around the last stable circular orbit. Based on this fact, we argue that 
the timing properties of the 7-ray signal could be determined by the processes in the direct vicinity 
of the supermassive black hole. 
Keywords: Galaxies: active 



1. INTRODUCTION 

Blazars are the subclass of radio-loud AGNs, for which 
the jet is closely aligned with the line of sight. They emit 
most of their power in the X- to 7-ray domain. This high- 
energy emission most probably originates in the jets and 
is produced via the Inverse Compton scattering of low- 
energy photons. 

Jets of the radio-loud AGNs are known to span over 
many orders of magnitude in distance, from the Astro- 
nomical Unit scale of the AGN central engine, the super- 
massive black hole (SMBH), up to the parsec to mega- 
parsec scales. The exact location of the "blazar emission 
zone" in the jet is not known. If the 7-ray flux origi- 
nates from within the Broad Line Region (BLR), where 
the radiation fields are intense, interactions of the ^-rays 
with low energy background photons leads to the e -pair 
produc tion thus prevent ing the 7-rays from escaping the 
source ([Liu fc Baill2QQ6[ ). This problem could be avoided 
if the site of production of the 7-rays is located outside 
from the BLR, parsecs away from the central SMBH. 

However, the possibility of location of the emission 
region at large distances poses a problem for the ex- 
planation of the observed variability timescales. Over 
the last few years, a number of blazars were found 
to exhibit variabili ty on very short timescales, down 
to se v eral minutes ([A haronian et al. 2007; Albert et al.l 
20071: lAbramows ki et al. 120101: JF oschin i et al.l 120111 : 
Neronov fc Vovkl 120111 : iSbarrato efaH 12011( 1. Simple 



causality arguments indicate that the short variability 
timescales could be best explained by the compactness 
of the 7-ray emission regions, which appears natural if 
the blazar emission zone is situated close to the cen- 
tral engine. In fact, 7-ray emission produced at the 
distances D ~ r 2 i?cE, where r ~ 3 — 30 is the bulk 
Lorentz factor of the jet and Rqe is the size of the 
blazar central engine, could be variable at the timescale 
as short as the light crossing time of the central engin e, 
tcE = Rce/c ~ 10 3 [Rce/2 AU] s (jCelotti et alJll998h . 
Although fast variability provides a strong argu- 
ment in favor of the 7-ray production within the 



BLR regions, close to the central engine, alterna- 
tive possibilities, usually invoking very compact emis- 
sion regions embedded in the larger scale jet have 
been cons idered to explain the fast variability phe- 
nomenon (Ghisellini & Tavecchio 2008; Ghisellini et al. 
120091 : iGiannios et al.l 120091 : iNaravanfc Piranll2012[ ). 

Accumulation of a larger sample of fast variable blazars 
is needed to distinguish between the two possibilities. If 
the fast variability is due to the location of the 7-ray 
emission region close to the blazar central engine, one 
generically expects that the shortest variability timescale 
is determined by the central engine's size. In this case 
the minimal variability timescales of blazars are expected 
to scale proportionally to the mass of the central SMBH. 
Otherwise, if the shortest timescales are related to the 
compact emission regions at large distances, formed as a 
result of the development of instabilities intrinsic to the 
jet, the minimal variability timescales are not expected 
to depend on the size of the SMBHs powering the jets. 

In this work we focus on the investigation of the min- 
imal variability timescales of a large number of blazars 
observed by the L arge Area Telescope (LAT) on board 
of Fermi satellite ([Atwood et al.l |2QQ9). LAT is operat- 
ing in the all-sky survey mode and is well suited for the 
monitoring of activity of blazars at different timescales 
in the 0.1-100 GeV energy band. At the same time, the 
relatively small effective area of LAT limits the photon 
statistics at the shortest minute-hour-day timescales ex- 
pected if the fastest variability scale is determined by 
the SMBH. Below we study the variability properties of 
a large number of blazars which are detected by LAT in 

\ 1 I K 

a two-year exposure (N olan et al.M2012[ ) and for which 
estimates of the mass of the SMBH exist. We demon- 
strate that although the minimal detected variability 
timescales are generically limited by the LAT sensitivity, 
for some ~ 10 brightest sources LAT sensitivity is suf- 
ficient to probe the fundamental SMBH timescale. For 
these brightest sources, variability at the timescale com- 
parable to the period of rotation around the last stable 
orbit is in fact observed. For 3 sources among these, 
the detected minimal variability timescales are found to 



be truly minimal, intrinsic to the sources. At the same 
time, we find that LAT sensitivity is currently not suffi- 
cient to confirm or rule out the correlation of the mini- 
mal variability timescales with the SMBH mass. We ar- 
gue that accumulation of larger exposure with LAT will 
increase the statistics of the bright blazar flares during 
which the sources exhibit the short-time-scale variability. 
This should open the possibility for testing the minimal 
variability timescale - black hole mass correlation. 

2. THE BLAZAR SAMPLE 

The Fermi Large Area Telescope (Fermi/LAT) is a 
space borne detector, sensitive to the 7 rays approx- 
imately from 100 MeV to 300 GeV. This instrument 
operates in all-sky scanning mode, completing the sur- 
vey in two revolutions (~ 3.2 hours). The Fe rmi/LAT 
two-year source catalogue (jNolan et al.l l2012f ) contains 



~ 1800 sources, among which ~ 1000 are associated with 
AGNs. 

For the purpose of this work, we select blazars from 
the LAT two-year catalogue for which the estimates of 
the black hole masses exist. The mass estimates mainly 
come from observations in the optical domain. Several 
alternative methods for the mass estimates are used. One 
possibility is to use the observed correlation between the 
SMBH mass and velocity dispersion in the central stel- 
lar bulge observed in the nearby galaxies. Measuring 
the bulge velocity dispersion through the spectroscopic 
observations of the blazar host g alaxy one could deduc e 
an estimate of the SMBH mass (jTremaine et al.ll2002f ). 
Otherwise, the black hole mass could be estimated from 
a combination of measurement of the width of the atomic 
line emission from the BLR and the time delays between 
the continu um and line emissi on (reverberation mapping 
technique) (jKaspi et al.ll2000f ). 
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1ES 0033+595 bzb 8.964 59.854 0.086 

1ES 0120+340 bzb 20.665 34.420 0.272 

Uy 0138-097 bzb 25.396 -9.481 1.034 

4C +15.05 agu 31.256 15.235 0.405 

PKS 0208-512 agu 32.696 -51.035 0.999 

3C 66A bzb 35.662 43.036 0.444 

4C +28.07 bzq 39.472 28.778 1.207 

AO 0235+164 bzb 39.675 16.624 0.940 

1H 0323+022 bzb 51.545 2.413 0.147 

PKS 0336-019 bzq 54.871 -1.744 0.852 

1H 0414+009 bzb 64.209 1.092 0.287 

PKS 0420-01 bzq 65.807 -1.342 0.915 

PKS 0440-00 bzq 70.686 -0.293 0.844 

PKS 0454-234 bzq 74.268 -23.427 1.003 

PKS 0454-46 bzq 74.041 -46.218 0.853 

4C -02.19 bzq 75.319 -1.931 2.286 

PKS 0521-36 agn 80.769 -36.469 0.055 

PKS 0528+134 bzq 82.714 13.553 2.070 

Uy 0537-286 bzq 84.846 -28.697 3.104 

PKS 0537-441 bzb 84.706 -44.084 0.894 

1ES 0647+250 bzb 102.699 25.096 0.203 

EXO 0706.1+5913 bzb 107.629 59.139 0.125 

S5 0716+714 bzb 110.476 71.350 0.300 

PKS 0735+17 bzb 114.524 17.703 0.424 

PKS 0754+100 bzb 119.288 9.963 0.266 

1ES 0806+524 bzb 122.460 52.308 0.137 

S4 0814+425 bzb 124.573 42.398 1.073 

PKS 0823+033 bzb 126.478 3.144 0.506 

B2 0827+24 bzq 127.646 24.122 0.940 

PKS 0829+046 bzb 127.987 4.486 0.174 

4C +71.07 bzq 130.420 70.880 2.172 

OJ 287 bzb 133.713 20.098 0.306 

S4 0917+44 bzq 140.236 44.697 2.189 

4C +55.17 bzq 149.433 55.382 0.895 

S4 0954+658 bzb 149.652 65.557 0.368 

1H 1013+498 bzb 153.788 49.432 0.212 

1ES 1028+511 bzb 157.761 50.895 0.360 

MKN 421 bzb 166.120 38.213 0.030 

1ES 1106+244 bzb 167.349 24.240 0.460 

PKS 1127-145 bzq 172.580 -14.803 1.187 

MKN 180 bzb 174.178 70.164 0.045 

4C +29.45 bzq 179.878 29.247 0.724 

ON 231 bzb 185.374 28.239 0.102 

4C +21.35 bzq 186.226 21.380 0.432 

3C 273 bzq 187.277 2.042 0.158 

PG 1246+586 bzb 192.063 58.350 0.847 

3C 279 bzq 194.042 -5.794 0.536 

Uy 1418+546 bzb 215.067 54.376 0.153 

PG 1424+240 bzb 216.760 23.795 0.160 

PG 1437+398 bzb 219.802 39.536 0.344 

1ES 1440+122 bzb 220.688 11.996 0.162 

BZB J1501+2238 bzb 225.275 22.639 0.235 

PKS 1510-089 bzq 228.207 -9.103 0.360 

1H 1515+660 bzb 229.517 65.437 0.702 

AP Lib bzb 229.435 -24.360 0.049 
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RGB J1534+372 


bzb 


233.861 


37.335 


0.143 


6.42 


8.13 + 0.70 
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57 


PG 1553+113 


bzb 


238.942 


11.190 


0.360 


— 


7.25 + 0.70 
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58 


PKS 1604+159 


bzb 


241.770 


15.876 


0.357 


— 


8.25 + 0.70 
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59 


4C +10.45 


bzq 


242.148 


10.494 


1.226 


— 


8.07 + 0.70 
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60 


OS 319 


bzq 


243.372 


34.159 


1.397 


— 


8.22 + 0.70 


3 


61 


4C +38.41 


bzq 


248.809 


38.171 


1.814 


4.81 


7.54 ± 0.70 
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62 


MKN 501 


bzb 


253.481 


39.763 


0.034 


— 


8.72 + 0.70 
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63 


NRAO 530 


bzq 


263.279 


-13.128 


0.902 


4.90 


8.07 + 0.70 


3 


64 


4C +51.37 


bzq 


265.060 


52.201 


1.375 


5.70 


7.97 + 0.70 


3 


65 


S4 1738+476 


bzq 


265.089 


47.645 


0.316 


— 


8.22 + 0.70 
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66 


1ES 1741+196 


bzb 


266.044 


19.579 


0.083 


— 


8.93 + 0.70 


1 


67 


Uy 1749+701 


bzb 


267.217 


70.111 


0.770 


4.68 


9.90 + 0.70 


1 


68 


OT 081 


bzb 


267.876 


9.640 


0.322 


5.47 


8.64 + 0.70 


1 


69 


S5 1803+784 


bzb 


270.147 


78.483 


0.680 


4.95 


8.82 + 0.70 


1 


70 


B2 1811+31 


bzb 


273.395 


31.722 


0.117 


6.26 


8.82 + 0.70 


1 


71 


4C +56.27 


bzb 


276.001 


56.838 


0.663 


6.60 


9.10 + 0.70 
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72 


PKS 1830-210 


bzq 


278.413 


-21.075 


2.507 


4.44 


8.05 + 0.70 


3 


73 


PKS 1933-400 


bzq 


294.320 


-39.933 


0.966 


— 


8.06 + 0.70 


3 


74 


1ES 1959+650 


bzb 


300.020 


65.157 


0.047 


— 


8.14 + 0.70 


1 


75 


Uy 2005-489 


bzb 


302.375 


-48.834 


0.071 


6.48 


8.57 + 0.14 
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76 


PKS 2052-47 


bzq 


314.068 


-47.255 


1.489 


5.57 


7.88 ± 0.70 


3 


77 


4C -02.81 


bzb 


323.466 


-1.906 


1.285 


5.88 


9.96 + 0.70 


1 


78 


S3 2141+17 


bzq 


325.879 


17.720 


0.211 


5.45 


8.14 + 0.34 
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79 


PKS 2149+173 


bzb 


328.102 


17.590 


0.000 


— 


7.78 ± 0.70 
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80 


PKS 2155-304 


bzb 


329.715 


-30.219 


0.117 


6.69 


7.60 ± 0.70 


3 


81 


BL Lac 


bzb 


330.707 


42.268 


0.069 


4.96 


8.58 + 0.70 


1 


82 


PKS 2201+04 


bzb 


331.164 


4.705 


0.027 


— 


7.76 + 0.13 
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83 


B3 2247+381 


bzb 


342.512 


38.419 


0.119 


6.15 


8.64 + 0.70 
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84 


3C 454.3 


bzq 


343.497 


16.153 


0.859 


2.94 


6.76 + 0.70 


3 


85 


BZB J2322+3436 


bzb 


350.672 


34.588 


0.098 


— 


6.76 + 0.70 


3 


86 


1ES 2344+514 


bzb 


356.759 


51.705 


0.044 


— 


8.57 + 0.70 


1 



Table 1 List of the sources used in this work. Bold marks sources found to 
be variable on the timescales comparable with the light-crossing time of the 

central black hole and /or the p e riod o f the last stable orbit. Sources 

classi fication is that of Nolan et al. (2012). Referenc es: (l)[W u et al. (2009), 

(2) IFalomo eTaLl (1 2003). (3) Fan et al. (2009), (4) Xie et al. (2005), (5) 



Paltani & Tiirler (2005), (6) Wagner (2008). 



In the case of blazars, which form the majority of the 
AGNs observed in the GeV domain, it is often the case 
that the velocity dispersion of the bulge can not be mea- 
sured due to the strong contamination with the emission 
of the AGN itself. In this case less direct estimates are 
used to infer the mass of the SMBH. These estimates are 
based on the so-called "Fundamental plane", that gives 
relation between the effective radius, average surface 
brightnes s and the velocity dispersion of the early-type 
galaxies (Fa lomo et al.ll2003l ). A similar correlation was 
found between the optical R-band lu minosity of the bulge 
and t he mass of the central SMBH (|McLure fc Dunlopl 
l2QQ2h . 

For the purpose of this wor k we c o mbined from the 
litera ture (|Falomo et all I2QQ2L 120031 : iPaltani fc Tiirlerl 
[20051 : iXie et al.ll2005HWu et al.ll2009h the list of 195 AGN 
with the SMBH masses measured using different meth- 
ods. We further selected those sources, that were ob- 
served by the Fermi/LAT telescope. Among those, only 
the sources without bright companions within 2°-radius 
around their positions were chosen for the analysis. Such 
selection criterion assures that the variability properties 
of the source of interest are not affected by the variability 
of nearby sources (the Point- Spread- Funciton of LAT is 
several degrees wide in the energy range below ~ 1 GeV 
Q). The resulting list of the 86 AGNs, used in present 
work, is given in table HJ 

3. DATA ANALYSIS 



1 http://www.slac.stanford.edu/exp/glast/groups/ 
/canda/lat_Performance.htm 



3.1. Light curves 

For the analysis we used the Fermi Science Tools 
v9r23pl^\ - a standard software package, provided by the 
Fermi collaboration to reduce the data, obtained by the 
Fermi/LAT. We also made use of the Pass 7 version of 
the photon files, which is stated by the Fermi collabo- 
ration as the best dataset for the Fermi/LAT analysifl 
We selected the class 2 photons and applied the cut on 
the zenith angle of 105°. 

To produce the light curves we used the gtbin tool, and 
computed exposures, corresponding to each bin, with the 
help of gtexposure. We fixed the spectral indices of the 
sources at the values, presented in the Fe rmi/LAT two- 
year source catalogue (jNolan et al.ll2012[ ). The energy 
range in our analyses was set to 100 MeV - 300 GeV. 
As we were not interested in the absolute normalization 
of flux, but were rather kin to detect variability, we set 
the software to collect photons within a small, 1° circle. 
This is much smaller than the 68% photon containment 
radius for 100 MeV, and is roughly equal to the same 
containment radius at 1 GeV. This way the analysis uses 
only the central, core part of the PSF, omitting its wings, 
thus minimizing contribution from the nearby sources. 
The light curves we built with the signal-to-noise ratio 
(SNR) binning, with the value of SNR in each bin fixed at 
4. Reducing the noise when the flux is low, such a binning 
allows to resolve finest details in the light curve during 
the flares and, thus, to probe the fastest variability of 

2 http://fermi.gsfc.nasa.gov/ssc/data/analysis/software/ 

3 http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/ 
/Pass7_usage.html 



the source. 

The above-mentioned procedure is equivalent to the 
usual aperture photometry with two exceptions. First is 
that we do not collect all the photons from the source, as 
we select only the narrow central core of the PSF for our 
analysis. The second difference comes from the fact, that 
we do not perform the background subtraction. As the 
background is not expected to be variable, this should 
not introduce any additional time scales into our results, 
but might limit the sensitivity to the fast variability of 
weak sources. 

3.2. Variability search 

To search for th e variability we utilize d the structure 
function approach (jSimonetti et al.lfl985h . The structure 
function (SF) of the light curve x(t) with the time lag r 
is defined as: 



SF x (T) = ((x{t + T)-x{t)f) 



(1) 



and represents the amount of variability at a given 
timescale r. It was first introduced in astronomy by 
iSimonetti et al.l (jl985f ) , and afterwards succ essfully used 
to me asure the variability of AGNs (e.g. iFavre et"aT1 
(|2QQ5h ). The advantage of the SF is that it works in 
the time domain, making it less sensitive to windowing 
problems than the Fourier analysis. 

There are a few important properties of the structure 
function, that are crucial in the frames of this paper. 
First of all, for very small values of r, far below the 
shortest variability timescale, the SF is dominated by the 
noise in the signal and, as the noise is independent from 
any timescale, the structure function here is constant. 
For the larger r, but still below the minimal variability 
timescale, the SF represents the small-scale linear trends 
in the signal and grows as r 2 . For the timescales, where 
the variability occurs, the change of the SF is defined 
by the type of the variability itself (for example, if a 
periodicity is present in the signal, the SF will show a 
dip at the corresponding period). Finally, for the largest 
values of r, larger than any variability timescale present 
in the signal, the structure function forms a plateau at 
the value, corresponding to the variance of the incoming 
signal. Thus the overall look of the structure function for 
a particular light curve is defined by the combination of 
the minimal and maximal variability timescales, the level 
of noise and the particular binning, used to produce the 
light curve. A deeper discussion of the SF properties and 
caveats was given by lEmmanoulopoulos et al.l (|2010f ) . 

The structure function also provides a direct way of 
measuring the amplitude of the variability at a given 
timescale. Ratio of values of SF at two different 
timescales is the square of the ratio of the corresponding 
amplitudes (see Eq. [I]). Probably the most convenient 
way of quantifying the variability amplitude is to com- 
pare the corresponding value of SF with its plateau value 
in the limit of large r, which represents the overall vari- 
ance of the light curve. 

To find the variability one has to compare the ob- 
served inclination slope of the logarithmically scaled 
SF (a=d\og(SF)/d\og(r)) with 2 (linear grows) and 
(below the minimal or above the maximal variability 
timescale). The linear growth, which not considered as 
variability here, may occur because of two different rea- 
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Figure 1. Example of the complete analysis cycle - 4C +21.35. 
Top left panel: the SNR-binned light curve. Top right panel: struc- 
ture function before (green) and after (red) the noise removal. 
Bottom panel: structure function logarithmic inclination index. 
Light blue lines represent the distribution of simulated SFs (central 
panel) and their indices (bottom panel), used to compute the error 
bars (95% confidence range). Minimal present variability timescale 
is 0.1 days. 

sons. First, expansion of almost any signal into the Tay- 
lor series leads to the dominant linear term for sufficiently 
small values of r, so generally one expects the slope of SF 
to be 2 for these values of r for any light curve. If, one 
the other hand, the light curve itself is a simple linear in- 
crease of flux, then there is no timescale to be associated 
with this growth, and the minimal variability timescale 
can not be derived. The apparent slope of the SF is af- 
fected by the noise, introduced due to the limited number 
of photons in each bin of the light curve. Before checking 
the slope of the SF one has to remove the contribution 
from this noise. In the this work we assumed that the 
noise is purely Poissonian, corresponding to the amount 
of counts in each time bin. To calculate its influence 
on the resulting SF, we produced a number of simulated 
light curves, for which the time bins were taken as in the 
initial one, and number of photons in each bin was taken 
from the Poisson distribution with the mean equal to the 
initial number of photons in this bin. This way we ex- 
pected to add to the light curve the very same noise as the 
one already present there, and, as a result, to double the 
value of noise plateau of the SF. The resulting difference 
between the averaged over the simulations SF and the 
initial one should be equal to value, added to the SF by 
the noise. Having subtracted this value from the initial 
SF, we obtained the noise-subtracted structure functions 
for each source. At the same time, the distribution of 
the values of simulated SFs for every given r defined the 
uncertainties of the structure function, derived from the 
original light curve. 

Uncertainties of the logarithmic inclination indices a of 
the structure function were also derived from these sim- 
ulations. Because the computed error bars of SF are cor- 
related (as every simulation resembles the original light 
curve) , we computed the inclination indices for each sim- 
ulated SF and then used their distribution for each value 
of r to compute the 95% confidence error bars for the 



indices of the initial structure function. 

We then compared all the indices of the SF with 2 and 
0, as it was described above. The value of r of the first 
bin, that, with the account of the error bars, would have 
an index different from these two reference values, was 
considered to be an estimate of the minimal variability 
timescale, present in the light curve. 

As an example, the light curve, structure function and 
structure function indices of 4C +21.35, obtained in this 
way, are shown in Fig. [TJ 

The minimal variability timescale, derived in this man- 
ner, may depend on the particular way the bins were 
chosen. Binning with a constant SNR essentially collects 
the photons in groups of N = SNR 2 , which is 16 in our 
case (SNR=4). Thus, in principle, we have a freedom 
of choosing these bins in 16 different ways. To account 
for this, for each source we produced 16 different light 
curves, each starting from the second photon of the pre- 
vious one, in this way going over all the possible combi- 
nations. The timescales, derived for these light curves, 
then gave us the uncertainty to which the minimal vari- 
ability is determined. If for any of those light curves 
variability was not found, we assumed that the source 
was not variable over the time interval covered by our 
analysis. 

4. RESULTS 

The derived minimal variability timescales, plotted 
versus the corresponding SMBH mass, are shown in 
Fig. [21 As it can be seen from this figure, the ma- 
jority of blazars in our sample vary on timescales sig- 
nificantly larger than the light-crossing times of their 
SMBHs (shown by the solid straight lines) 
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where R g = GMbh/c 2 and a = Jbh/Mbhc 2 is defined 
by the angular momentum Jbh of the black hole and lies 
in the range < a < R g . 

A small fraction (some ~ 10%) of the sources demon- 
strate variability timescales as small the period of the 
last stable orbit (shown by the dashed lines in Fig. [2]) 
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From Fig. [2] one can notice that different blazar types, 
BL Lacs and Flat Spectrum Radio Quasars (FSRQ), ex- 
hibit variabilities at different timescales, with FSRQs be- 
ing systematically variable at shorter timescales. It is not 
clear a-priori, if this reflects the difference in the mech- 
anisms of 7-ray emission in these types of objects or is 
related to a selection bias or an instrumental effect in 
LAT. 

It is also not clear a-priori if the sources which do not 
reveal variability on the timescales as short as the black 
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Figure 2. Measured minimal variability timescales versus the 
mass of the central black hole for the AGNs from table [T] Solid 
lines represent the light-crossing times for maximally rotating (red) 
and non-rotating (blue) black holes, dashed lines - the same for the 
period of the last stable prograde orbit. 

hole light crossing time and/or the period of rotation 
around the last stable orbit are just not variable on these 
timescales or the sensitivity of the LAT is not sufficient 
for the detection of variability on such short timescales. 

In order to investigate the sensitivity limit of LAT, we 
plot in Fig. [3] the dependence of the minimal variability 
timescale on the peak count rate in the LAT. From this 
figure one could see a trend that the minimal measured 
variability timescale decreases with the increasing peak 
count rate. This could be readily explained by the fact 
that detection of variability at a given timescale requires 
sufficient signal statistics. The minimal statistics in our 
case is 16 photons in one time bin. The minimal de- 
tectable variability time could not be shorter than the 
time necessary for collection of 32 photons (there should 
be at least two time bins to compare the fluxes). This 
minimal detectable scale is shown by the dashed line 
in Fig. 03 One could see that the detected variability 
timescales follow the minimal detectable timescale line. 
Thus, the non-detection of variability of the signal from 
most of the blazars on the SMBH timescales is not re- 
lated to the absence of such variability, but is rather due 
to the sensitivity limit of LAT. 

From the same figure one could see that the difference 
in the variability timescales of FSRQs and BL Lacs is due 
to the fact that FSRQs are on average brighter than BL 
Lacs in the LAT energy band. Higher count rates of FS- 
RQs allow the detection of shorter variability timescales 
for these sources. 

One can also notice from Fig. [3j that for 3 blazars, 
namely 3C 273, 4C +21.35 and PKS 1510-089, the sen- 
sitivity of Fermi/LAT allows the detection of an order of 
magnitude faster variability, than really observed. This 
suggests that the found minimal variability timescales of 
these sources correspond to the truly minimal timescales, 
defined by the physics of their emission regions. Analysis 
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Figure 3. Measured minimal variability timescales versus the 
maximal count rate found in present analysis. Dashed line rep- 
resents the minimal variability timescale measurable with a given 
count rate. Sources found to have variability on the timescales 
comparable to/shorter than the period of the last stable orbit are 
marked with bold. 

results for these sources are briefly summarized below. 

3C 273 was variable by factor ~ 100 in terms of 
flux over the period of the Fermi/LAT observations. 
This highly variable source is also one of the bright- 
est in our sample, thus allowing to study the variabil- 
ity at the smallest timescales. This is demonstrated 
in Fig. [4j where one of the bright flares, happened on 
the 21th of September 2009, is shown. As it can be 
seen from this figure, brightness of this flare ensures 
that the sufficient amount of photons per bin is achieved 
even for several hours long intervals, so variability down 
to a few hours could be, in principle, detectable for 
this particular light curve, if the variability amplitude 
is large enough, whereas our analysis indicated, that 
only ~ day timescales are present (Fig. [3j). This can 
be seen from Fig. [5j which shows the structure func- 
tion for this particular light curve. Minimal variability 
timescale, found by our analysis of the entire light curve, 
is 0.58lo'.45 days. The corresponding light-crossing time 
of the central black hole and period of rotation at the last 
stable orbit can be estimated from equations [2] and G3 
ti c ~ 0-3 -T- 0.6 days and P(r m i n ) ~ 1.7 -r- 14 days, de- 
pending on the black hole rotation momentum. 

A curious feature, seen in Fig. [5j is a dip at the 
timescale of ~ 5 days, which may indicate (quasi)periodic 
changes in the light curve. Indeed, this behaviour is 
seen in the light curve, with the most prominent episode 
shown in Fig. [6] If being associated with the emission of 
some blobs of gas, just about to fall onto the black hole, 
this gives us the period of rotation on the last stable or- 
bit around the central SMBH. As the exact value of this 
period depends on the black hole's rotation momentum, 
similar measurements (with a better statistics) could be 
used to measure the latter. 

4C +21.35 is another highly variable source with 
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Figure 4. A part of the light curve of 3C 273 around the 21th of 
September 2009 (JD 2455096), with the bright flare in it. 




Figure 5. Structure function of the light curve of 3C 273, par- 
tially shown in Fig. [6] Error bars correspond to 68% confidence 
ranges. The structure functions from the simulated light curves, 
used to determine the error bars, are shown with light blue. Mini- 
mal detected variability timescale - « 1.18 days - is marked with 
the vertical dashed line. A dip in the structure function is seen at 
r ~ 5 days. 

bright flares (see Fig. [1]). The minimal variability 
timescale, determined by our analysis, is 0.050_]q'q24 days 

or 1.20^0.58 nr - Visual inspection of the light curve sug- 
gests that it took place during the brightest flare on the 
30th of April 2010. As Fermi/LAT operates in the all-sky 
scanning mode with the scanning period of 3.2 hr, this 
measurements of the variability timescales below 3 hr 
may be affected by the fact that telescope did not ob- 
serve the source constantly. Consequently, it would be 
more conservative to say that the minimal detected vari- 
ability timescale lies between 0.5 hr, which is the average 
duration of the observation of a given patch of the sky, 
and 3.2 hr. 

PKS 1510-089 has produced a number of bright flares 
during the period of the Fermi/LAT operation, with the 
strongest one taking place on the 19th of October 2011. 
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Figure 6. A part of the light curve of 3C 273 in September- 
November 2009, which shows a possible qausi-periodicity in it. 

The minimal variability timescale, derived through the 
SF analysis, is 0.111q!o6 days or 2-6-i'.5 hr. Similarly to 
the case of 4C +21.35, we visually inspected the light 
curve and concluded that this estimate corresponds to 
the above mentioned flare. As the derived time scale 
is smaller than 3.2 hr, this estimate is subject to the 
same uncertainty, as the estimate for 4C +21.35. At r ~ 
50 days structure function exhibits a dip, that continues 
up to the timescales comparable with the duration of 
the light curve, which is the largest timescale we can test 
with our analysis. 

5. DISCUSSION 

Minimal variability timescale provides information on 
the size and location of the 7-ray emission region and/or 
on the physical process responsible for the emission. In 
the case of blazars, the high-energy 7-ray emission might 
originate from a "blazar emission zone" in the jet at large 
distance from the central SMBH or, alternatively, from 
the direct vicinity of the black hole. We have proposed a 
method to distinguish between these two scenaria, via a 
study of correlation of the minimal variability timescale 
of 7-ray emission with the mass of the SMBH. In the 
former case no such correlation is expected, because the 
shortest variability timescales are determined by the size 
scales and/or Lorentz factors of small inhomogeneities 
in the blazar jets (possibly formed via development of 
plasma instabilities or magnetic reconnection events). In 
the latter scenario, a correlation is generically expected 
because the shortest variability timescale is determined 
by the size of the blazar central engine. 

We have attempted to verify the existence of the " min- 
imal timescale" - "black hole mass" correlation, using 
the data of Fermi/LAT. Our analysis shows that the 
LAT sensitivity is sufficient for detection of the variabil- 
ity timescales as short as the SMBH light crossing time 
and/or period of rotation around the last stable circular 
orbit for a number of bright flaring blazars, see Fig. [2j 
From this figure one could see that some ~ 10 blazars 
(out of 86 considered in the analysis) exhibit variability 
on the "fundamental" SMBH timescales. At the same 
time, none of the sources is variable on the timescales 
shorter than the black hole light crossing time. 



From the same figure one could see that most of the 
sources are variable at least on the timescales within an 
order of magnitude from the period of rotation around 
the last stable orbit around Schwarzschild black hole. 
For most of the sources the minimal variability timescale 
measured by LAT is determined by the low statistics of 
the signal in LAT, rather than by the physics of the 7- 
ray emission intrinsic to the source. This is clear from 
Fig. [3j in which a trend of the decrease of the minimal 
measured variability timescale with the increase of the 
peak flux is observed. Only three sources in our sample 
(3C 273, 4C +21.35 and PKS 1510-089) deviate from the 
formal sensitivity limit due to the finite 7-ray countrate 
in LAT, shown by the dashed line in Fig. 03 

In these three sources the measured minimal variabil- 
ity timescale might be determined by the physics of the 
source, rather than by the LAT performance. In this case 
the measurement of the minimal variability timescale 
would provide an important information on the nature of 
the 7-ray emission. Indeed, the estimates of the masses 
of black holes in 4C +21.35, PKS 1510-089 and 3C 273 
differ by more than an order of magnitude, but in all the 
three sources the minimal variability timescale derived 
from the LAT data is close to the light crossing time 
of the SMBH. This provides a support for the hypothe- 
sis of relation of the minimal variability timescale to the 
size of the blazar central engine. It is clear, however, 
that accumulation of larger sample of blazars, for which 
the intrinsic minimal variability timescale is measured, 
is needed before definitive conclusions on the validity of 
this hypothesis could be drawn. 

One possible caveat of the measurement of the minimal 
variability timescale in 4C +21.35, PKS 1510-089 and 3C 
273 could be in the peculiarity of the observation mode 
of LAT. LAT images each particular direction on the sky 
for the time spans of ~ 0.5 hr once every 3.2 hours. De- 
tection of variability on the timescales between 0.5 and 
3.2 hours could, therefore, be affected by the irregularity 
of exposure by LAT. In the cases of PKS 1510-089 and 
4C +21.35, the minimal detectable variability timescales 
are shorter than 0.5 hr so that, in principle, the vari- 
ability could be measured within a single 0.5 hr expo- 
sure (with several time bins contained in this exposure). 
However, in the case of 3C 273, the minimal measurable 
variability timescale is close to 0.5 hr, so that the vari- 
ability on this timescale might be missed if a particular 
flare of the source is not fully contained within a single 
0.5 hr exposure. At the same time, 3C 273 produced 
a number of bright flares over the three year exposure 
with LAT, so that it is perhaps unlikely that half- an- hour 
scale variability has been missed in all the flares. Clari- 
fication of the issue of the instrumental 3.2 hr and 0.5 hr 
scales for the search of the minimal variability timescales 
will become possible with accumulation of more flaring 
sources with minimal detectable variability timescales in 
this range. 

If the observed minimal variability timescales are in- 
deed related to the timescales of the SMBH in the blazar 
central engine, they provide a direct measurement of 
the size of the central engine, without an uncertainty 
related to the Dopp ler factor of the emitting plasma 
(|Neronov et al.ll200ct ). Indeed, as the black hole and the 
observer are located in the same rest frame, the variabil- 
ity produced by the black hole will be seen by observer 



at the very same timescale, independently of what hap- 
pens on the way with the medium, carrying the emission. 
Thus, the only difference between the intrinsic and ob- 
served timescales comes from the cosmological redshift 
of the source, which enlarges the observed timescales by 
(1 + z) with respect to the intrinsic ones. Thus, 7-ray 
observations could potentially provide independent es- 
timates of the masses of SMBHs in the cores of AGN, 
providing a important tool in the context of the study 
of the formation and evolution of the SMBHs. Indeed, 
in our analysis only 86 sources with known estimates of 
the SMBH masses were included (see T able 1). However, 
LAT AGN catalog ([Nolan et all 12012( 1 includes a much 
larger number of variable sources, for which the measure- 
ments of the minimal variability timescales could be ob- 
tained, resulting in the estimates (or upper limits) on the 
masses of the SMBHs. The 7-ray tool for the estimate of 
the black hole masses would be usable only after accu- 
mulation of larger statistics of the measurements of the 
intrinsic minimal variability timescales of 7-ray emission 
from blazars and verification of the possible correlation 
between the black hole mass and variability timescale 
shown in Fig. O 

As it was mentioned before, current instruments, such 
as Fermi/LAT, allow only for the brightest objects to be 
studied at the sufficient level of details. The question 
whether the minimal variability timescale is indeed lim- 
ited by the size of central black hole, thus, might need 
to await for the new generation of 7-ray instruments to 
become operational. 
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